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Abstract. The planetary nebulae He 2-436 and Wray 16- 
423 in the Sagittarius dwarf galaxy appear to result from 
nearly twin stars, except that third-dredge-up carbon is 
more abundant in He 2-436. A thorough photoionization- 
model analysis implies that ratios Ne/0, S/0 and Ar/0 
are significantly smaller in He 2-436, indicative of third- 
dredge-up oxygen enrichment. The enrichment of oxygen 
with respect to carbon is (7±4)%. Excess nitrogen in Wray 
16-423 suggests third dredge-up of late-CN-cycle products 
even in these low-mass, intermediate-metallicity stars. 

Key words: Nuclear reactions, nucleosynthesis, abun- 
dances - (ISM:) planetary nebulae: general - (ISM:) plan- 
etary nebulae: individual: He 2-436 and Wray 16-423 - 
Galaxies: dwarf - Galaxies: individual: Sagittarius 



1. Introduction 



It is usually assumed {e.g. Henry 1989) that the oxygen 



abundance of a planetary nebula (PN) reflects the metal- 
licity Z of the interstellar medium (ISM) in which the 
parent star formed, inasmuch as the conversion of ^^O 
into ^"^N due to advanced CNO processing is unimpor- 
tant, a condition fulfilled in low-mass stars {M < 2 Mq). 
Low-mass stars bring to their surface nitrogen (through 
first dredge- up), produced by CN-cycle conversion of car- 
bon, and carbon (third dredge-up), produced by 3-a fu- 
sion along the A symptotic Giant Branch (AGB, e.g. Iben 
& Renzini 1983| ). During He-burning, ^^N is transformed 
notably into ^^Ne and some ^^C into ^^O. Thus third 
dredge-up is likely to bring freshly synthesized oxygen to 
the surface of low- m ass stars, as suggested by Boothroyd 
& Sackmann (1988, here BS88) and, in a preliminary form. 



by Herwig et al. (200C, here HBD), who invoked intershell 
dredge- up ("fourth dredge- up", Iben 199£) by diffusive 
overshoot. Exposed deep layers of hydrogen-deficient post- 
AGB stars appear to have quite similar concentrations of 



oxygen and carbon (Koesterke & Hamann 1997; Herwig 
et al. 1999) but the amount of fresh oxygen expelled from 
AGB stars is unknown. The oxygen abundance is difficult 
to determine in the a tmosphere of evolved stars {e.g. Ful- 
bright & Kraft 1999) and shows large scatt er in PNe of a 
given galaxy {e.g. Leisy & Dennefeld 1996 , here LD96). 

In order to reveal the oxygen enrichment of PNe, one 
must seek for circumstances enabling the detection of a 
small excess of oxygen. Two PNe, He 2-436 and Wray 16- 
423, w ere fo und (Zijlstra & Walsh [1996| , here ZW96; Walsh 
et al. 1997) to belong to the Sagittarius dwarf galaxy 
(Ibata et al. 1995). Spatial and kinematic properties point 
to a common origin for these PNe. Photoionization models 



(Dudziak et al. 2000, here DPZW) indicate nearly iden- 
tical depletions with respect to solar for all elements be- 
yond nitrogen [(—0.55 ± 0.07) dex]. Both PN nuclei are 
early-type [WC] stars belonging to nearly the same (H- 
burning) evolutionar y trac k M=(1.2 ± 0.1)Mq, Z=0.004 
(Vassiliadis & Wood 1994), consistent with their precur- 
sors being the intermediate-age carbon stars of Sagittar- 
ius. As expected, third dredge-up yield was large in these 
intermediate-^ stars (Marigo et al. 1999), but even larger 
in He 2-436. If both stars were born in the same star forma- 
tion episode with identical abundances, advantage can be 
taken of their low initial abundances and different third- 
dredge-up yields to study the influence of this process 
upon oxygen and other elements by comparison of the 
respective abundances. Here, unlike the general study of 
DPZW, emphasis is on abundance ratios and differential 
effects. Pequignot et al. ( 2000| ) gave a preliminary account. 
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2. The Sagittarius dwarf galaxy PN models 

Spectra of Wray 16-423 and He 2-436 were secured in iden- 
tical conditions. Their small apparent size led to global 
spectra of the highest reproducibility and relevance for 
modeling. The many theoretical line ratios (H i, He i, for- 
bidden multiplets) allowed to check that statistical errors 
on measured fluxes were realistic, the few discrepancies 
being indeed similar in both PNe. Thus, in a direct com- 
parison, systematic errors arc likely to cancel out. 

In the photoionization model analysis of DPZW, the 
error bars attached to the abundances of these PNe were 
estimated from comparing two independent codes: they 
represented, at some "Icr level", the uncertainties left af- 
ter evaluating the origin of the differences between codes 
and trying to include systematic errors. Here, only one 
of the photoionization codes (code NEBU; e.g. Morisset 
& Pequignot 1996 ) is used, considering only the statisti- 
cal error on line fluxes. Because of the neglect of system- 
atic errors, the intervals are much reduced. Conversely, 
because of the "plasticity" allowed by the rather large 
number of free parameters and the thorough exploration 
of all possible solutions, the intervals can be enlarged. The 
abundance ratios, which were not considered in the error 
analysis of DPZW, are determined with greater accuracy 
than the individual abundances, as shown below. 

Current [WC]-atmosphere models are not reliable: in 
order not to miss solutions because of restrictive assump- 
tions, black-body stars with a discontinuity at the He+ 
limit were employed as a three-parameter sequence and an 
equivalence was established with the NLTE model stellar 
atmosphere (effective temperature Tcff), from the grid of 



Clegg & Middlemass (1987). Similar free-parameter sets. 



procedures and convergence criteria were adopted to ob- 
tain two-sector models for both PNe. The models matched 
essentially all line fluxes within Icr errors, that is no more 
than a few percent of line strength for most basic lines. 
The two discrepancies left, one of them involving [O i] and 
the other one [Ariv]/[Ariii], turned out to be identical for 
both PNe, suggesting systematic effects. 

It was possible to generate several models of each 
PN satisfying the very strong and numerous observa- 
tional constraints, as three key line fluxes (Cii 426.7, 
[On] 372.6 and [Sii] 406. 8nm) were not very accurately 
determined. Given a Tcff for the central star, the model 
parameters are constrained by imaging and/or radio con- 
tinuum data and by (accurate) optical line fluxes. A no- 
ticeable aspect is the global energy balance. The cool- 
ing rate, mainly determined by the unknown C ill] and 
Civ UV line fluxes (in addition to [Oiii]), depends on 
the carbon abundance, to which C ii is proportional. For 
Wray 16-423, the computed Cii flux is always in agree- 
ment with observation, stronger limitations being brought 
by other lines. For He 2-436, considering different models 
with increasing density in the dense sector (an incomplete 
inner shell), both the [O iii]500.7nm/[O iii]436.3nTO ratio 



and the [O iii]500.7nm flux tend to decrease, but can be 
reconciled with observation by increasing the oxygen and 
carbon abundances. As free parameters ensure sufficient 
flexibility, the uniqueness of the solution breaks down and 
the boundaries of the set of acceptable models (and abun- 
dances) are now determined by the C il flux. 

In Fig. 1 are shown for both PNe the ranges of al- 
lowed carbon and oxygen abundances versus the Toff's. 
Varying Teff makes the predicted [Sii] 406. Snm and [On] 
372. 6nm fluxes vary relative to other lines in models. For 
He 2-436, low or high Teff's are forbidden as the [Sii] or 
[O ii] fluxes become too large. Relaxing somewhat the er- 
ror bars attached to these lines increases the allowed range 
of Teff without much changing the range of acceptable 
abundances. As for Wray 16-423, all elements but carbon 
have abundances virtually independent of Teff. Low Teff's 
(low C abundances) are forbidden by two independent 
[Six] lines being predicted too weak; high Teff's by one 
accurately observed [S ii] line becoming too large. Higher 
Teff's are also excluded by the fact the He"*" discontinuity, 
constrained by He ii 468. 6nm, would be much larger than 
in atmosphere models. Also the shape of the model nebula 
would then deviate from observation. 

The oxygen abundances of the He 2-436 models en- 
compass the well defined value of Wray 16-423. Neverthe- 
less oxygen tends to be larger in He 2-436. As a check 
of this tendency, Wray 16-423 models were run with the 
oxygen abundance increased by 10%. No satisfactory so- 
lution could be found. Re-convergence on all parameters 
was performed until all lines except [O ill] 500.7nrn and 
[O III] 436. 3nm were again satisfactory and the discrep- 
ancy on both [O III] line fluxes was reduced to a mini- 
mum, typically 5a. Other elemental abundances increased 
by 9±1% (except C: 18%). An He 2-436 model of given 
Teff was then re-converged to a solution, though with the 
condition that the computed [Olll] line fluxes should re- 
produce the discrepancy met in the previous Wray 16-423 
model. In this new He 2-436 model, most abundances were 
shifted upwards relative to the standard solution by 8±1% 
(except N: 3% and C: 18%), close to the shift imposed on 
oxygen in the Wray 16-423 model. This numerical exper- 
iment suggests that systematic errors may substantially 
shift the derived abundances but not influence conclusions 
based on differences between the two PN models. 

Table 1 lists the absolute helium and oxygen abun- 
dances by number and other abundances relative to oxy- 
gen. For all elements but carbon, the quoted abundance 
"uncertainty" is obtained by adding linearly the full abun- 
dance spread of acceptable models to the la error of the 
best determined line (the models were constrained to re- 
produce this line exactly). Percent differences on abun- 
dance ratios and la error bars (Table 1, col. 4) clearly 
show that Ne/0, S/0 and Ar/0 are smaller in He 2-436 
than in Wray 16-423 (3.5, 1.8 and l.Aa confidence levels re- 
spectively, hence 4.2cr collectively). Given that sulfur and 
argon are not synthesized in low-mass stars and that neon 
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Table 2. Evolution of abundances* in stellar envelopes 
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Fig. 1. Acceptable abundances of C and O by number 
versus T^s of central stars according to model nebulae 
for He 2-436 and Wray 16-423. Best-model values are in- 
dicated. The shaded areas correspond to He 2-436 models 
having the same sulfur abundance as Wray 16-423 within 
uncertainties (Sect. 3). For Toff/IO^K = 6.5, lO'^xO/H in- 
creases from 24 to 28 as lO^xC/H increases from 104 to 
135. For Teff/IO^'K ^ 7.5, the O and C values are 24 and 
140 respectively. 

Table 1. Abundance ratios by number in Sagittarius PNe 



Ratio 


Wray 16-423 


He 2-436 


diff. (%) 


He/H 


0.108±0.002 


0.108±0.003 


0.0±3.6 


O/H 


(21.4±0.2)xl0"^ 


(24.0±4.0)xl0"^ 


-11.5±16. 


C/0 


3.22±0.23 


5.0±0.8 


-43.3±17. 


N/0 


0.213±0.015 


0.116±0.006 


59.0±8.7 


Ne/0 


0.168±0.003 


0.147±0.005 


13.3±3.8 


S/O 


0.0201±0.0008 


0.0162±0.0016 


24.0±13. 


Ar/0 


0.0042±0.0004 


0.0029±0.0007 


36.6±26. 



should hardly be affected either, this very significant trend 
must be due to a larger oxygen abundance in He 2-436. 
An explanation involving an abundance fluctuation in the 
ISM with, e.g., an excess specific to oxygen at the birth 
place of the parent star of He 2-436 is unlikely, as these 
stars are twins in most respects (ZW96, DPZW) and the 
Ne/0 ratio is remarkably stable, at ~ 0.18, in the ISM of 



even different galaxies {e.g. Henry 1989). Thus the most 
plausible interpretation of the excess is in terms of oxy- 
gen production, with the precursor star of He 2-436 having 
been more efficient than that of Wray 16-423. Since He2- 
436 shows a larger carbon excess - the signature of third 
dredge-up - it is inferred that this process was also able 
to bring freshly synthesized oxygen to the star surface. 



Stage: 


Init 


1st drge-up 


3rd drge-up 


3rd drge-up 


Elmt 




+ GBP" 


(Wray) 


(He 2) 


He 


8250 


9000 


10800±200 


10800±300 





(17) 


16.5±3.2 


21.4±0.2 


26.0±2.0 


C 


4 


2.5 


69.±5. 


122. ±18. 


N 


1 


2.5 


4.6±0.3 


3.0±0.3 


Ne 




3.4 


3.59±0.07 


3.8±0.3 


S 




0.43 


0.43±0.02 


0.43±0.02 


Ar 




0.08 


0.09±0.01 


0.07±0.02 



'Abundances by number in units H = 10"^ 



"Cool Bottom Processing (Boothroyd & Sackmann 199£) 



3. Discussion 

In Table 2 is sketched a possible evolution of the stellar 
envelope composition, assuming that (1) initial composi- 
tions were identical in the Wray 16-423 and He 2-436 stars, 
(2) sulfur abundance was constant, keeping the value de- 
termined for Wray 16-423 within uncertainties (col. 4 of 
Table 2), and (3) oxygen abundance varied linearly with 
that of carbon during third dredge- up (cols. 4, 5). Sulfur 
is preferred to argon as a reference constant-abundance el- 
ement, in view of the moderate scatter of the sulfur abun- 
dance determinations in models and the difficulty met in 
accounting for the argon ionization (DPZW). Scatter is 
less for neon, but its nucleogenic status may be less clear- 
cut and the variation of Ne/0 from Wray 16-423 to He2- 
436, although significant, is smaller. Reliable sulfur abun- 
dances are difficult to obtain in PNe, yet this approach 
is justified considering that models are based on homoge- 
neous procedures and used only in a differential manner. 

Once the sulfur abundance is specified for He 2-436, 
the range of possible oxygen abundances is narrowed: as- 
suming the errors are independent, S/O (Table 1, col. 3) 
and S (Table 2, col. 4) imply 23.6 < 10^ x (O/H) < 28, an 
interval that excludes the O/H value of Wray 16-423. Ow- 
ing to the linear increase of S with O in the model results, 
the minimum O/H is in fact slightly larger (Table 2, col. 5; 
Fig. 1, shaded areas). Let AO and AC be the increments 
by number of the oxygen and carbon abundances due to 
third dredge-up. Taking into account correlations between 
C and O in He 2-436 models (Fig. 1: the lowest C goes with 
the lowest O) and assumption (3) above, the difference in 
abundances between the two PNe shows that extreme val- 
ues for AO/AC are 2.4/76=0.03 and 6.8/61=0.11. 

The extrapolated main-sequence abundance of oxygen 
(Table 2, col. 3 or 2), obtained assuming carbon is ini- 
tially much less abundant than in the PNe, is used only to 
determine Z and relevant initial abundances (see below). 

The present result for AO/ AC seems to fit with old 
predictions of BS88 for flash-produced intershell 0/C in 
the region mixed after convection but this may be fortu- 
itous: th ese s tellar models are uncertain, as indicated by 
Herwig ( 200C ) who obtained models for intermediate- mass 
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high-Z stars using an exponential diffusive overshoot for- 
malism. Based on the same concept, HBD obtained pre- 
liminary models for low-Z stars, suggesting that surface 
0/C reaches asymptotically the intershell value, ~ 1/3. 
Thus, considering tendencies shown by recent star models, 
the present coarse estimate would seem to imply that rela- 
tively little oxygen should be third-dredged-up. Nonethe- 
less a PN is formed after many ejection episodes and its 
abundance may be far from reflecting asymptotic values. 

By contrast, neon and argon appear constant within 
uncertainties and the lower abundance of nitrogen in He 2- 
436 cannot be due to extensive conversion of this ele- 
ment into heavier ones. The initial (Table 2, col. 2) and 
first-dredge-up (col. 3) helium, carbon and nitrogen abun- 
dances are from Boothroyd & Sackmann ( |1999 ) assuming 
M=1.2 Mq and Z=0.004 (in keeping with extrapolated 
oxygen). Second dredge-up does not occur in low-mass 
stars, but these authors take into account deep-circulation 
mixing leading to enhanced ^'^C/^^C, as observed along 
the red giant branch, and to further conversion of initial 
carbon into nitrogen. This so-called "Cool Bottom Pro- 
cessing", sensitive to Z and M, can bring the nitrogen 
abundance quite close to that of He 2-436 but perhaps not 
to the larger one of Wray 16-423 (Table 2). Another source 
of nitrogen is required. 

Mixing of ^H with fresh ^^C occurs in low-mass AGB 
stars to produce ^^C, s- nuclei etc, but the hydrodynamics 
involved is not understood (Lattanzio & Frost 1997). It 
may not be excluded that some ^"^C can find its way up to 
the hydrogen shell or that insertion of H-rich pockets into 
the He/C-rich intershell is, on occasion, efficient enough 
to allow CN-cycle completion. In this way, low-mass AGB 
stars could synthesize some ^^N from fresh ^^C and mimic 
intermediate-mass stars in which the so-called "Hot Bot- 
tom Burning" , together with third dredge-up, is believed 
to lead to (N-rich) Type I PNe. Within the uncertain- 
ties of the definition (Torres-Peimbert & Peimbert 1997 ) , 
Wray 16-423 would in fact be marginally Type 1 in the 
Large Magellanic Cloud (LMC). The indication here is 
that the dredge-up of additional nitrogen was more effec- 
tive in Wray 16-423, wliere carbon is less overabundant. 

According to BS88, the He/C ratio in the carbon pock- 
ets that will be dredged-up levels off rapidly at He/C 8 
by number, compatible with the "observed" third-dredge- 
up AHe/AC, which lies anywhere between 15 and —15 
(Table 2, cols. 4, 5), but not with the large jump of he- 
lium between first dredge-up and Wray 16-423. This may 
suggest, among other possibilities, that the initial helium 
abundance adopted here for these stars was too small. 

Abundances are comparable in the PNe of Sagittarius 
and LMC (LD96). The abundance scatter in the LMC is 
likely due to the coexistence of several populations, but 
Ne/0 is expected to be relatively stable, except possibly 
for differential effects of the kind revealed here. ^From 
data of LD96, it is found that log(Ne/0) is -0.61±0.10 
and -0.83±0.10 for Type-1 and non-Type-1 PNe respec- 



tively. These determinations bracket the mean H il-region 
value, -0.74 (see LD96), and appear marginally different. 
Assuming that neon best traces the overall LMC evolu- 
tion, this suggests that oxygen is relatively depleted in 
Type 1, at least as a consequence of second dredge-up in 
intermediate-mass stars (see also Dopita & Meathering- 



ham, 1991), but enhanced in non-Type 1 by third dredge- 



up in accord with present finding in Sagittarius PNe. 

Abundances are lower in the SMC than in the LMC. 
Nitrogen-rich PNe of the SMC (type of LD96) are 
found among extremely low Z stars, presumably of low 
initial mass. Also some halo PNe have N/0 > 1 (Howard 



et al. 1997), at variance with the well-known association 



of Type I PNe with intermediate-mass stars, but maybe 
another manifestation of what is occuring in Sagittarius: 
^'^N is enhanced by third dredge-up in low-mass AGB stars 
and, due to a combination of efficiency and titration effect, 
N/0 can reach large values for sufficiently small Z. 

The discussion of the Sagittarius PN models (DPZW) 
makes it clear that improvement of the critical carbon 
abundance can be expected from more complete observa- 
tions (UV, far red). Nevertheless the present results can 
already provide some constraints on the new generation of 
models for the evolution of low-mass \ow-Z stars and lead 
to a re-interpretation of trends shown by PN abundances 
in low-Z environments. The production of oxygen by low- 
mass stars (~ 2xlO^"'Af0/star in Sagittarius) may have 
to be considered in galactic-evolution models. 
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